Intrinsic spin Hall effect at oxide interfaces: a simple model 
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"^ Abstract 

T-H An asymmetric triangular potential well provides the simplest model for the confinement of 

'~~^ mobile electrons at the interface between two insulating oxides, such as LaAlOa and SrTiOa 

(LAO/STO). These electrons have been recently shown to exhibit a large spin-orbit coupling of the 






^ Rashba type, i.e., linear in the in-plane momentum. In this paper we study the intrinsic spin Hall 

C effect due to Rashba coupling in an asymmetric triangular potential well. Besides splitting each 

C^ subband into two branches of opposite helicity, the spin-orbit interaction causes the wave function 

I in the direction perpendicular to the plane of the quantum well (the z direction) to depend on 

^ the in-plane wave vector k. In contrast to the extreme asymmetric case, i.e., the wedge-shaped 

I quantum well, for which the intrinsic spin Hall effect vanishes due to vertex corrections, we find 

^""1 that the asymmetric well supports a non- vanishing intrinsic spin Hall conductivity, which increases 
> 

^^ in magnitude as the well becomes more symmetric. The spin Hall conductivity is found to be pro- 






portional to the square of the spin-orbit coupling constant and, in the limit of low carrier density. 



■^ depends only on the effective mass renormalization associated with the k-dependence of the wave 

O 

^^ functions in the z direction. Its origin lies in the non- vanishing matrix elements of the spin current 

^ between subbands corresponding to different states of quantized motion perpendicular to the plane 

/\ of the well. 

J3 



I. INTRODUCTION 

The spin Hall effeclP has been a topic of great interest in the last decade^"^, and it has now 
become a mainstream technique for the manipulation of spins in spintronic devices^^"^. Its 
signature is the appearance of a current of z-spin in the ^/-direction following the application 
of an electric field in the x direction.'^^^^The inverse effect, i.e., the generation of a transverse 
electric field by an injected spin current has also been observed.'^SHffil 

It is by now clear that the spin Hall effect results from an intricate competition of several 
mechanisms.'^ In all cases one needs a "sink of momentum" , usually provided by impurities 
or phonons, in order to attain a steady-state response to the applied electric field. The 
spin Hall effect in a crystalline solid can be described as extrinsic or intrinsic depending 
on whether it is driven by spin-orbit interaction with impurities (extrinsic case) or with the 
atomic cores of the regular lattice (intrinsic case). Early studies focused on the analytically 
solvable model of a two-dimensional electron gas in a wedge-shaped quantum well with 
Rashba spin-orbit coupling.'^Si This model is relevant to (001) GaAs quantum wells in the 
presence of an electric field perpendicular to the plane of the electrons. In this model, 
the intrinsic and extrinsic components of the effect are distinguished by symmetry - the 
former being even under a reversal of the sign of the spin-orbit coupling constant, while 
the latter is odd. It was soon realized that the intrinsic spin Hall conductivity for this 
model vanishes exactly,'23H35] ^^^ ^j^g extrinsic spin Hall conductivity is suppressed when the 
Dyakonov-Perel spin relaxation rate becomes comparable to or exceeds the Elliott- Yafet 
spin relaxation rate.'^^"^ On the other hand, the vanishing of the spin Hall conductivity 
has been recognized to be a peculiarity of this model.'^ Intrinsic spin Hall conductivity has 
been predicted and observed in two-dimensional hole gases, and, most remarkably, in three 
dimensional centro-symmetric d-band metals (e.g. Pt).'^ 

Recently, a high mobility two-dimensional electron gas (2DEG) of tunable density 
has been observed at the interface between two insulating oxides, such as LaAlOs and 
SrTiOa.'^^"^ A large spin-orbit splitting of the Rashba form, i.e., AE(k) = hak, where k is 
the wave vector in the interfacial plane and ha can be as large as 5 x 10^^ eV A, has been 
observed in the proximity of a superconducting transition.HS In view of this large spin-orbit 
splitting, the interfacial 2DEG seems an excellent candidate for the observation of the spin 
Hall effect. However, one must be wary of the fact that vertex corrections tend to suppress 



the spin Hall conductivity of systems with linear-in-A; spin-orbit splitting. 

In this paper, we introduce a very simple, analytically solvable model, which we hope 
will help clarify the essential ingredients of the intrinsic spin Hall effect at oxide interfaces. 
The model consists of an asymmetric triangular quantum well in which we have a spin-orbit 
interaction of the Rashba type, acting only on the right half [z > 0) of the well. Thus, the 
Hamiltonian is 

where Q{z) is the Heaviside step function {Q{z) = 1 for ^ > and Q{z) = otherwise); p = 
{px,Py) is the momentum in the plane of the quantum well; z is the coordinate perpendicular 
to the plane, and Ac is the effective "Compton wavelength" , which controls the strength of the 
spin-orbit coupling in the relevant conduction band of the quantum well. In semiconductors 
like GaAs Ac is known to scale inversely to the cube of the fundamental band gap and 
amounts to a few Angstroms. In oxide materials it is somewhat smaller: Ac — 0.7 A as can 
be inferred from the observed value of ha = X^eE ~ 5 x 10^^ eV A assuming an electric 
field of the order of 1 V/A. The model potential is 

i Fz z>0 
V{z) = , (2) 

I -rFz z <0 

where r > 1 is our asymmetry parameter and F = eE, E being the magnitude of the electric 
field in the z direction and e the absolute value of the electron charge. 

In spite of its simplicity, this model gives a reasonable description of electrons bound at 
the interface of two insulating oxides, such as SrTiOs and LaAlOa.'^^Sll We have in mind an 
n-type interface, which is equivalent to a sheet of positive charge at 2; = 0. According to 
band-structure calculations^ the electrons that neutralize this sheet of positive charge reside 
primarily on the SrTiOs side {z > 0) where both the band gap and the electric field are 
smaller. On the LaAlOs side {z < 0) the electric field is larger, due to reduced electrostatic 
screening, and this is the effect we try to capture with the parameter r > 1. Further, the 
spin-orbit interaction within the conduction band is largely determined by the spin-orbit 
interaction of the "B" ion within the perovskite formula ABO3; the spin-orbit interaction of 
the Al orbitals is negligible compared with that of the Ti orbitals, so the Rashba spin-orbit 
coupling will be much smaller on the LaAlOs side than on the SrTiOa side. This is the 



situation modeled with the Heaviside function; the spin-orbit couphng is significant only for 
2> 0. 

The eigenfunctions of the hamiltonian ([I]) have the form 

gikr 1 / 1 \ 

where A is the area of the interface, k = {k^, ky) is the in-plane wave vector, r is the position 
in the interfacial plane and z is the coordinate perpendicular to the plane, ^k is the angle 
between k and the x axis. These states are classified by a subband index n = 0, 1, 2.., which 
plays the role of principal quantum number, an in-plane wave vector k, and a helicity index, 
A = +1 or —1 which determines the form of the spin-dependent part of the wave function. 
The most interesting feature of this model is that the subband wave functions fnkxiz) 
depend on k. This feature is crucial to the existence of a non-vanishing intrinsic spin Hall 
conductivity. This can be seen most clearly by applying to the present model the standard 
argument for the vanishing of the intrinsic spin Hall conductivity in the Rashba model.SSl 
According to Eq. ([I| the time derivative of ay is 

^y = {[H, <yy] = -2^V\z)e{z)pya, . (4) 

The expectation value of a time derivative must vanish in a steady state, hence 

{V'iz)Q{z)pya,) = Y, Nnkx{nkX\V'iz)Q{z)kya,\nkX) = , (5) 

nkA 

where A^nkA is the average occupation numbers of jrikA) in the given non-equilibrium state. 
If at this point we were allowed to factor the average into a product of {V'{z)Q{z)) and 
(pyCTz), we could immediately conclude that the spin Hall current, being proportional to 
{Py<^z), is zero. This argument works in the Rashba model because the z-dependent part of 
the wave function does not depend on k or A. In the present case, however, the z-dependent 
wave functions fnkxiz) do depend on k and A, creating a correlation between V'{z)Q{z) and 
PyCTz- Then, we can no longer assert that the spin Hall current vanishes. In fact, writing 

{V'{z)e{z)pyaz) = {V\z)e{z)){pyaz) + (A[K'(^)e(z)]A[p,a,]) (6) 

where A [A] represents the fluctuation of a quantity relative to its average, we can conclude 
that 

(A[n^)e(^)]Ab^ 



Our calculations confirm this. The shape of the confining potential in our model is 
controlled by the asymmetry parameter r. For r — )■ oo the electrons are entirely confined to 
the right half of the quantum well. In this limit we recover the Rashba model. The subband 
wave functions become independent of k, because the spin-orbit coupling is independent of 
z in the region of space in which the electrons move. The spin Hall conductivity vanishes, 
when the "vertex correction" is duly taken into account. 

For finite r additional contributions of order a^ appear from the /c-dependence of the 
transverse wave functions. Here a is the dimensionless spin-orbit coupling constant 

«-f^V> (8) 



where 

1/3 



(—) 

\2mFj 



(9) 



is the natural length scale of our model. Since i is of the order of a few Angstroms for oxide 
interfaces, while Ac — lA, we see that a is somewhat smaller than 1, but not orders of mag- 
nitude smaller. We believe that a multi-band effect, proportional to a^, is also responsible 
for the intrinsic spin Hall conductivity predicted and observed in centro-symmetric metals 
like Pt.'^ Indeed, we have verified that the a^ spin Hall conductivity is present in our model 
even if we make the spin-orbit coupling symmetric and set r = 1, thus simulating a centro- 
symmetric material. In this case the subbands become doubly degenerate with respect to 
the helicity index, but the a^ contribution to the spin Hall conductivity is still present. 

Assuming, as discussed above, that our model gives a reasonably good description of 
electrons at oxide interfaces, the results of our work imply that the intrinsic spin Hall effect 
will be an effect of order a^ (as opposed to a°) and be crucially dependent on the energy 
separation between the transverse subbands. 

This paper is organized as follows. In Section II we discuss the analytic solution of the 
model. In Section HI we calculate the inter-subband contributions to the SHC as a function 
of r. In Section IV we calculate the intra-subband contribution to the spin Hall effect - first 
without including vertex corrections (where it is found to be quite large and independent 
of a ) and then with vertex corrections (the correct approach). In the latter case, only a 
contribution proportional to a^ survives. 



II. SOLUTION OF THE MODEL 



We express all quantities in units derived from the natural length scale £ defined in 
Eq.(N), i.e., lengths in units of £, momenta in units of h/d ~ lA~^, energies in units of 




Ft = 2^;^ ~ 1 eV. All quantities in the following treatment are therefore dimensionless, 
unless noted otherwise. The dimensionless hamiltonian takes the form 

^ = 5Z {^' - ^ + ^('^) - ae{z)v\z)[Kay - kyaA (10) 

where k = p/h in physical units, and k = p in the present reduced units. Here we have 
defined 

viz) = I (11) 

I —rz z < . 

We note that kx, ky and kxCy — kyax are compatible constants of the motion, the latter 

with eigenvalues Xk with A = ±1 and eigenstates of the form 

(12) 

where 6'k is the angle between k and the x axis. Therefore we classify the eigenstates by 
quantum numbers k (in-plane wave vector) and A = ±1 (helicity). The wave functions for 
the n-th subband (n = 0, 1, 2... in order of increasing energy) have the form given in Eq. ([s]) 
where fnkx{z) is the solution of the Schrodinger equation 

- fnkxiz) + [v{z) - XaQ{z)kv'{z)]fnkxiz) = enkxfnkx{z) (13) 

with 

i I z>0, 
v'{z) = (14) 

I -r z <0. 

The eigenvalues e„kA are arranged in increasing order, 

cofcA < eifcA < e2fcA < ••• (15) 

The complete energy is 

Enkx = k'^ + e„kA • (16) 



The solution of the Schrodinger equation (13) for generic n is 



I -^;(;-^'-^--^f) (z > 0) 
M^) = zy<{ ..(^;g^^/3/ - \, , (17) 



Ai{-zr'^/^-e^:^r-^/^) /^ ^ qn 




FIG. 1: Wave functions for the three lowest subbands plotted vs z for ak = 0.2 for r = 2. Notice 
the difference between A = 1 and A = — 1 wavefunctions with the same n, denoted by dashed (red) 
and solid (blue) lines respectively. Notice that the vertical axis has been shifted to separate clearly 
the wave functions with different n. In reality, all wave functions tend to for \z\ — )• oo. 



where A%{x) is the Airy function and Z is the normaUzation constant. Notice that this wave 
function is continuous at 2 = 0, with /kA(0) = Z. The energy ekA is determined by imposing 
the continuity of the derivative at z = 0: 



A^{-tkx - Aafc) _ _ ^i^ Ai{-tkxr-^l^) 



A%{-ek\ - Xak) 



Az(-efcAr-2/3) 



(18) 



where the Ai'(z) = 



dAi(z) 
dz 



For small k the subband energies and wave functions can be obtained from standard 



perturbation theory. Let us denote by fn{,z) the solutions of Eq. (13) at A; = and by e„ 



the corresponding energies (clearly, these are independent of a and A). Then, we have 



fnkx{z) = fn{z) + Xak ^ ^^^^^fn>{z) + ... , 



where 



and 



Pn'n = {fn'\v'{z)Q{z)\fn) = I fn'{z)fn{z)dz, 

'o 



CnfcA = CnO " Xakpnn + [akf ^ "^ h .. 



n'^n 



(19) 

(20) 

(21) 



n '-n' 




Elk- 




Eok+ 



FIG. 2: The two lowest subband energies e^kA plotted as functions of ak for r = 2. Notice the 
linear splitting between opposite helicities. The middle line in each case is the average of the 
two helicity bands and shows the (barely perceptible) curvature associated with the effective mass 
renormalization. 

Bands of equal n and opposite helicities cross at A; = 0. The intra-subband splitting is 
linear in k and proportional to Pnn- The quadratic correction indicates the emergence of a 
spin-orbit-induced effective mass m*, which is related to the bare band mass m by 



m 
in* 






n'^n 



(22) 



All of the above features are confirmed by numerical calculations (see figure below). 
The following facts are emphasized, which will play a key role in the calculation of the 
spin Hall conductivity: 

1. The subband wave functions are functions of the product ak. Wave functions of oppo- 
site helicities fnk+{z) and fnk-{z) are generally different (see Fig. [I]). This difference 
disappears only for ak = Q (i.e., either /c = or a = 0). This situation is in stark 
contrast with the conventional Rashba model, in which the subband wave functions 
are independent of a, k or A. 

2. Subband wave functions of the same k and opposite helicities are not mutually orthog- 



onal, even when their subband indices are different. From Eq. (19) we see that 

Pn'n 



{fnk\\fn'k-\) = -2akX- 



^n ^n/ 



n ^ n . 



(23) 



up to terms of order (ak)"^. On the other hand, for n = n' we have 



ifnkxlfnkx) — 1 
{fnkx\fnk-x) = l-2iakyYl 



Pn' 



^n ^n' 



(24) 



To derive these overlaps it is essential to take into account second-order corrections 
to fnkx, ensuring that the latter remains normalized to 1 to second order in ak. More 
generally, we have 



ifnkxlfnk'x) = I - ^Hk - k')]' J2 
{fnkx\fnk'-x) = 1 - l[a{k + k')]' Y. 



n'^n 



^n ^n' 



Pn' 



Cn Cri/ 



(25) 



These formulas will be needed later in the analysis of vertex corrections. 



III. INTRINSIC SPIN HALL EFFECT 



2 

yx 



In this section we begin the calculation of the intrinsic spin Hall conductivity [ctsh] 
which connects the jy = kya^ component of the spin current to an electric field applied in 
the X direction. The real part of the spin Hall conductivity is given by the formula 



\'^SH_ 



yx 



e 
A 



E 



Qm 



dip 



nkA 



dA' 



dil), 



nkA 



dAl 



(26) 



nkA, occ 

where the sum runs over the occupied states. We will assume in what follows that only states 
in the lowest subband (i.e., the subband with ra = 0) are significantly populated in the ground 
state. We use perturbation theory to calculate the variation of the wave functions ipnkx to 
the application of infinitesimal vector potentials that couple to k^ and ky respectively, the 
first shifting k^ — )■ fc^ + ^°, the second shifting ky ^- ky + Ayffz (the hamiltonian must be 
properly symmetrized after the second replacement, to preserve hermiticity). It is evident 
that these homogeneous perturbations do not mix wave functions of different k. For the 
asymmetric well, it is not necessary to use degenerate perturbation theory. Instead, the 
usual non-degenerate theory suffices. The zeroth order wavefunctions are given by Eq. dsj). 



where the subband wave functions fnkx{z) are given by Eq. (17), evaluated at the appropriate 



values of the energy. In order to apply our formula ( 26 ) we need to find the single-particle 



wave functions to first order in the applied potentials A^^^Ay. We note that the application 
of y4° modifies the wavefunction at a given k, A in the following manner: 



^nkA[A°](r) 




fniix{z)e 



jkr 



(27) 



where k = k + yl°.x is the shifted wave vector. Thus, specializing to the ri = subband, we 
find 



f^^okA(r) 



dAl 



X 



A0=0 



sin6'k 

V2k 







/okA(^)e^'^-'- + 



V2 \ iXe^s^ I dk 



Next, we want to calculate 



dtl>ok\{r) 



dAl 



Ag=0 



To do this, we observe that the first order change 



in the wavefunction due to the coupling to Ay is 

^V'okAlr) = 2 Alky Y, ^^""''^°''^ ^„.kA(r) 
^ eokA - e„kA 



(29) 



where A is a short-hand for —A. The reason why only states of helicity —A appear in the 
variation of ■i/'okA is that the Ay field couples to the spin current operator kya^, which fiips 
the helicity index. Writing ky as k sin ^k yields 



c^V'okA(r) 



dA' 



2k sin 6'k / 



Ag=0 



(/nkAl/okA) 
CokA — e„kA 



ij. 



nkAl 



(30) 



Lastly, we combine Eqs. (28) and (30) for the derivatives of the wave function into the 



general formula (26) for the spin Hall conductivity. Only the first term on the right hand 



side of Eq. (28) contributes, due to the vanishing of the angular average of sin 6'k cos ^k- 



Taking into account the fact that the lowest subband states are populated in the ground 
state up to fc = kp+ for helicity +1 and up to /c = kp- for helicity —1 we arrive, after some 
simple algebra, to the following expression. 



[^SH_ 



yx 



e 



E 



L^o 



eOk+ — Cnk- Jo ^Ok- ~ enk+ 



(31) 



All the quantities that appear in the square brackets of this equation are dimensionless. 
We now examine separately the inter-subband contributions {n ^ 0) and the intra-subband 
contributions [n = 0) to the spin Hall conductivity. 
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FIG. 3: Intersubband contribution to the spin Hall conductivity in units of 



ea^kp 



(see Eq. ([32j)). 

Data points have been connected to enhance visibility of the trend. Contributions up to n = 5 
have been included. Notice the logarithmic scale for r. 



Inter-subband contribution 



Let us first focus on the contribution of the terms with n ^ 0. It is safe to assume 



that kp^ and kp^ are sufficiently small to justify the use of the small-fc expression (23) for 



the overlap between wave functions of opposite helicity. In this approximation we can also 
ignore the difference between kp^ and kp-, i.e., we set kp^ = kp^ = kp. Further, we can 
ignore the /c-dependence of the subband energies, i.e., we set e„kA — £«• Then, making use 



of Eq. (23) and carrying out the k integral, we arrive at: 

1 7,4 V^ \Pn-o 



a 



inter] 2 
SH \yx 



-a 



IT 



n^O 



(eo 



(32) 



All the quantities on the right hand side, except e, are dimensionless. Thus kp = \/27mP, 
where n is the two-dimensional electronic density and i is the length scale of the confining 
potential. For ^ = 3 A and n = 10^^ cm~^ we have kp ~ 0.75. Although derived for 
the clean limit - i.e., for frequencies that are much larger than the inverse of the electron 
impurity scattering time 1/r (while still much smaller than the Fermi energy Ep) - this 
result is essentially exact for the d.c. transport regime as well. Vertex corrections do not 
modify it significantly, since they are of the order (Ar)^^ ^ 1 where r is the momentum 
relaxation time and A is the energy separation between the subbands. 

The numerically calculated inter-subband contribution to the spin Hall conductivity is 
plotted (in units of -^ ) in Fig. Ill A We note that it is present for all values of r, but it 



11 



vanishes in the Rashba hmit (r — )■ oo) because pno vanishes in that hmit. We also notice that 
this contribution to the spin Hall conductivity vanishes in the limit of low carrier density, 
i.e., for kp -^ 0. 

B. Intra-subband contribution 



Let us now consider the intra-subband contribution 

e f'^^ \{f-o\Uo)? 



a 



intral z 



47r ./fc^_ eok+ - eok- 



rCt drC . 



(33) 



'SH \yx - 

where we have exploited the fact that |(/ok-|/ok+)P = |(/ok+|/ok-)P to combine the two 



integrals on the right hand side of Eq. (31) into a single one 



The wave vectors kp^ and kp- are determined by the equations 



Kp: -\- Kp_ ZK 



where kp = 27ini'^, and 



Eokp++ — EQkp_ 



F 5 



We expand the energy up to third order in a: 

Eokx = k^ + ei{ak)X + e2{akf + e3(aA;)^A 
where 



ei 



-Poo 

n^O 



bo '^ 



JOn\ 



1/772 



E 



POnPnmPmO 



(en -eo)(er 



eoj 



Poo 



E 



\p 



On I 



eoj 



_n,m^O n^O 

Then the solution for the Fermi wave vectors can be written in the form 



where 



kpx = kp + kiaX + k2(y^ + k^a^^X. 



ki 
k2 
k. 



ei 

2 




8kp 




6162 


esk]. 


2 


2 
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(34) 
(35) 

(36) 

(37) 
(38) 

(39) 

(40) 

(41) 
(42) 
(43) 



Armed with these results we easily find that 



and 



kp^ — kp- — —aei — a kpC^ + a 6162 



eok+ - eok- - 2akei + 2{akfe:>, . 



We also know that (see Eq. (24)) 



\{fou+\hu-)?^l~^{ak)^Y. 



PnO 



Co — Cn 



Combining these results we arrive at 



cr 



„intra] z 
SH \yx 



e m 



in m 



l-AiakpfYl 



PnO 



Co ~ Cn 



(44) 



(45) 



(46) 



(47) 



where we have used the fact that 2 



m ^^ rn_ 
m* m 



. This result is striking, because its leading 



order is independent of both a and kp. However, as we now show, this result is completely 
invalidated by vertex corrections, which cancel the a° term and leave us with an a^ term, 
which remains finite when kp ^ 0. 

C. Vertex corrections 



In this subsection we discuss the vertex corrections to the spin Hall conductivity. To this 
end, it is practical to resort to the standard diagrammatic formulation in terms of the Green 
functions. The zero-temperature Kubo formula for the spin Hall conductivity reads 



Wsh];^ = -e\un- r |^Tr{%^G(e+)7.G(e_)}, 



(48) 



where the trace is taken in the basis of the exact eigenstates |/) = ^Jnkxijiz). The energies 
in the Green functions are e± = e ± uj/2, and e > is the unit charge. In the reduced 
units, the current vertices read 7^ = kya^ and ^x = — a6(z)cr^. Notice that v'{z) = 1 for 
2; > has been absorbed into a and only the anomalous spin-dependent part of ^x has been 
considered, since the "regular" part does not contribute to the spin Hall conductivity. In 
the basis of the exact eigenstates the Green function is diagonal and reads 

1 



G,(e) 



e-Ei + iO+sgn{Ei 



(49) 
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El = Enkx = k"^ + e„kA being the energy eigenvalues. By inserting the resolution of the 
identity, I = ^^ 1 0(^1 5 twice under the trace and performing the integral over e, we get 

By exploiting the fact that i{l\'yy\l'){l'\'yx\l) is real, after using the Kramers-Kronig relations, 
we can rewrite the spin Hall conductivity in the form 

miy\n{inx\i) — . _ ,2 — • (51) 

By making the shifts k^ ^ k^ + A^ and ky ^ ky + AyCr^ we obtain the bare current and 
spin-current vertices: 

7j-^|| -*.-". (52) 

where the ~ in the first equation reminds us that the "regular" part of the charge current 
vertex has been omitted. 

The relevant matrix elements of these current vertices are the ones between states of 
opposite helicity: 

(nkA|7^|n'kA) = kyifn^xlL'-kx) (53) 

(nkA|7^|n'kA) ~ -i^^lKkA - ^„'kA)(/nkA|/n'kA)- (54) 



Eq. ( 53 ) can be derived almost by inspection. Eq. ( 54 ) can be obtained by using the eigenvalue 
equation for the functions /nkA(-2)- 

Following the procedure described in Ref. [Ml we calculate the renormalized current vertex 
Fa; according to the equations 

% = -ae{z)cry + ^-i- J2 G^' ^K Gt' ■ (55) 

Here A^o = l/47r is the dimensionless density of states in the absence of spin orbit interaction. 
The superscripts R and A stand for retarded and advanced Green functions, respectively. 
The first equation represents the ladder resummation for an effective vertex 73^., which is 

14 



defined by the second equation. Notice that the first term in the expression for 7^; is the 



bare vertex of Eq. (54). In the purely two-dimensional Rashba model with no ^-dependence, 
one sees that the second term on the right hand side of the second equation cancels the first. 
This is the famous vertex cancellation. To see this explicitly one must project the above 
equation into the spin states |A), the projection over the plane wave states already being 
done. 

To extend the treatment to the present case, the projection must be made over the states 
\nX). Within the approximation of disorder with no 2;-dependence of the impurity potential. 



the vertex equations are not changed. The second of the Eqs. (55) becomes 



(nkA|7^|n'kA') = ~{nkX\ae{z)ay\n'kX) 

+ ^^ E {nkX\n,k'X^)Gl^,,^ 2^ Gi^,,^{mk'X,\n'kX'). (56) 

nik'Ai 

The matrix elements (rakA|r2ik'Ai) and {riik' Xi\n'kX') are those of the impurity potential. 
Explicitly we have 

(nkAlmk'Ai) = ^(/„kA|/n,k'A,)(l + AAie'(^><'-^'')) (57) 

(mk'Ailn'kA') = ^(/„,k'Aj/n'kA') (1 + A'Aie-^^^'^'-^'')) . (58) 

By observing that k'^ = k' cos 6'^', one can perform the integration over the direction of k' in 



Eq.(56) 



1 r ^ (1 + AAie^^^'-'-'''')) cos(^k') (1 + X'X^e-'^'-'-'-A = ^ (Ae-^'' + X'e''-) . (59) 
4 Jq zn 8 



We now rewrite Eq. ( 56 ) for the matrix elements necessary for the evaluation of the spin Hall 



conductivity. By using Eq.(54) 



(OkA|7.|nkA) = -i^^(i?0kA-i5„kA)(/0kA|/„kA) 

~ i g^jY ^ 2^ -^l(/okA|/nik'Ai)(/nik'Ai|/„kA)'^nik'Ai 2/i;' G^^k'Ai- (60) 
nik'Ai 

For weak disorder, we may take the limit r — )■ oo and perform the integral over k' to get 
(OkA|7.|nkA) = -i^^(i?okA-i?„kA)(/okA|/„kA) 

~ i — ^T7 7,-^1 (/okA I fmkFm Ai Ai ) {fnikp„^ x^ Ai I /nkx) ^Fn-^Xi A^mAi, (61) 



2iVo ^ 

niA 
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kpniXi and A^mAi being the Fermi moinentuin and the density of states at the chemical 
potential for the rii-th subband with helicity Ai. Let us consider first the intra-band matrix 



elements, n = 0, so that Eq.(|61|) becomes 

(0kA|7,|0kA) = 

/ , ^l{fokx\fnikp„^XT^l)ifnikF„,x,\i\fo]ix) ^PniXi iVniAi, (62) 



i (-^OkA - -E'okA)(/okA|/okA) 



— 1- 



2k 
.sin(6'k)A 



2N<, 



niAi 



Now we show how this equation can be evaluated in a small a expansion up to the third 
order. We begin by observing that by performing the sum over Ai, always one of the overlap 
factors yields a Kronecker 5„j o because of the orthonormality of the wave functions /„kA- 



Furthermore, because of Eqs.(25), we can replace kp^x^ with k with an accuracy a^. In this 
way the overlap factor between states of the same helicity can be approximated with unity, 
while the other is common to the bare part of the vertex (the first term on the right hand 



side of Eq.(62)). Hence we are reduced to evaluating 

.sin(6'k)A 



(0kA|7.|0kA) 
ifokxlfokx) 



-l^^(i?OkA 



^Okx) 



2Nn 



7 ^ Ai kpn^Xi Nn-iXi- 



(63) 



The densities of states at the Fermi level in the lowest subband are evaluated from the 
formula 

(64) 



JVa ^ 2fcA 

^0 IVfc-EofcAl/cA 

where kx is the solution of the equation Eq^x = Ep. By using the equations (36 39) for the 



energy and (40 43) for the Fermi wave vectors, we obt 


ain 






No 2kp 


'6162 Seskp el 
Akp 2 16fc|_ 




The quantity relevant for the vertex correction is 




kxNx kp (^ X ei 2 3 

— 7^ = — s 1 — Aa- ha- — 

2No 2 [ kp kp 


el 62 ^F 
Skp 2 


kp 


'ei62 
2 


e^ky 
2 


The vertex correction is then 




kxNx k-xN-x _ 


\fyp. -^' 


^\n,^ 


I Pi Pn Pnh/ 


l^ 





(65) 



(66) 



2No 2Na 



while the bare vertex is 



2k 



{Eokx - Eokx) = Aa6i + Aa^63fc^ 



(67) 



(68) 
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where, in the last term, k can be replaced by kp at the required level of accuracy. By 
summing the above two expressions, the effective vertex reads 



(0kA|7x.|0kA) = -iAsin(^k)(/okA|/okA)«' (2eie2 - e^kl) 



(69) 



The above expression must replace Eq.(54) when evaluating the intra-band spin Hall con- 



ductivity and one can also safely neglect the ladder resummation in the first equation of (55) 
due to the r — > oo limit. Hence, to the accuracy we are working it is enough to multiply the 



bare bubble intra-band spin Hall conductivity of Eq.(47) by the factor — «^ (2eie2 — e^kj 



Notice that in obtaining Eq.(47), the momentum integral must be evaluated with accuracy 
a^ in order to keep the corrections up to order a^ in the SHC The presence of the corrected 
vertex, which is already of the order a^ smaller than the bare one, allows the evaluation of 
the integral with accuracy a. Notice also that in the Rashba limit, r — t- oo, the coefficients 
Pon vanish, and hence 62 and 63 both vanish, yielding the vertex cancellation as expected. 



The intra-band spin Hall conductivity of Eq. (|47|) must then be replaced by 

POnPnmPmO 



^intral ^ 
\yx 



\ intra~\ - 



-a 



■i-n 






br '^ 



JQn 



I |_ ^F 

eo Poo 



,n,m^O 



(cji ~ eojl^m ~ Co) 



n^O 



bo '^ 



JOn\ 



,en - eoj 



which is accurate up to terms of order a^. This can also be written as 



[ intral 



a'^kp 



yx 




E 



POnPnmPmO 



eo)(er 



eo 



E 



bonP 
n — Co) 



(70) 



(71) 



The numerically calculated intra-subband contribution to the spin Hall conductivity from 



Eq. (70) is plotted (in units of ^^- ) in Fig. 4. We note that it is present for all values of r, 
but it vanishes in the Rashba limit (r — )■ 00) because pnm vanishes in that limit. 

Vertex corrections, in principle, are present also for the inter-band matrix elements con- 
trolling the inter-band spin Hall conductivity. In the case of the matrix elements between 



the occupied and unoccupied sub-bands, with accuracy of order a, Eq.(61) becomes 

.sin(6'k)A 

niAi 



(OkA|7x.|nkA) = -isin(6'k)Aapon-i- 



2Na 



y , M{fokpx\fnikpXi){fnikpXi\fnkp\)kFniXiNniXi- 

(72) 



By performing the sum over Ai and using again the orthonormality of the wave functions 
fnkx, one obtains 

Poo + Pnn 



(OkA|7x|?2kA) = -isin(6'k)Aapon - isin(6'k)apo^ 



eo 



(73) 
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ia^ln 



FIG. 4: Intrasubband contribution to the spin Hall conductivity in units of ^^ for kp = l(see 
Eq. (|32|)). Data points have been connected to enhance visibility of the trend. Notice the logarith- 
mic scale for r. 



We see that the first vertex correction for the inter-band matrix elements contains as a 
denominator the energy separation. For large enough separation, this correction can be 
neglected. It is interesting to note that the energy separation in the denominator appears 
because of the overlap factors due to the impurity potential scattering. The physical ori- 
gin of these processes is the following. Upon scattering from an impurity an electron can 
make a transition to an unoccupied subband, where afterwards changes its spin direction by 
precessing in the Rashba field. Clearly these processes cost energy and yield a small correc- 
tion to the spin-dependent inter-band matrix elements of the current vertex. As a result, 



Eq.(32) for the inter-band spin Hall conductivity is not changed much, either qualitatively 



or quantitatively, by the vertex corrections. 



Eqs. (32) and (71) are the main results of this paper, showing a spin Hall conductivity 



of order a^. Notice that, in the low-density limit /ci? — )■ only the intra-band term survives, 
reducing to the simple form 

I--); - ~ (l - ^) . (74) 

which is entirely controlled by the spin-orbit-induced effective mass. Fig. 5 examines the 



relative importance of the two terms in the curly brackets of Eq. (71) at kp = 1- The 



effective mass term (first term) is plotted in red, the second term is plotted in black, and 
their sum is in blue (this coincides with the result plotted in Fig. 4). We see that the two 



terms in the curly brackets of Eq. (71) have opposite signs, but the first term dominates 
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?0.01 



-0.03 



FIG. 5: Comparison of the first and the second term in the curly brackets of Eq. (71) for the 



intra-band spin Hah conductivity at kp = 1. First term, denoted by I, is in red, second term, 
denoted by II, is in black; their sum is in blue. Notice that the two terms have opposite signs and 
the first term dominates. 

even at kp as large as 1. 



IV. CONCLUSION 

We have developed a simple model for the 2DEG that exists at the interface between two 
oxides. Neglecting band structure effects and including only an asymmetric wedge-shaped 
potential that binds the electrons to the interface, and the spin-orbit interaction associated 
with it, we have calculated the intrinsic spin Hall conductivity in the high- mobility limit. 
After a careful consideration of vertex corrections to the single-bubble result, we have found 
that the intrinsic SHC is finite and, in the low-carrier density limit, has the simple form 



of Eq. (74), which crucially depends on the mass renormalization associated with the k- 
dependence of the wave functions in the z direction. This effect, which is of order a^, 
vanishes in the standard Rashba model, which is the r — )■ oo limit of the present model. 
Finally, the numerical value of the SHC that we calculate is of the order of ^0.08a^. This 
should be compared with the values measured in the two-dimensional electron gas in GaAsjl^S 
which can be expressed as ^0.05a^. Thus our values are of a magnitude comparable with 
those reported in GaAs if a is of order 1 as reported in the literature.'^Sl 

Our model clearly does not include the spin-orbit interaction that is built into the two- 
dimensional band structure of the interfacial electrons, coming from the spin-orbit inter- 
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action with the atomic cores. The latter can be extracted from a tight-binding model 
calculation for SrTiOa, taking into account the fact that the interfacial electrons live almost 
entirely on the STO side of the LAO/STO interface.SSl From the form of the in-plane wave 
function in the relevant conduction band (arising from Ti d-orbitals split by crystal fields 
and spin-orbit interaction with the atomic core) one can extract an effective spin-orbit cou- 
pled hamiltonian, and the intrinsic spin Hall conductivity can be calculated. A comparison 
between the SHC calculated in this paper and that obtained from a tight-binding calculation 
of the spin-orbit coupled band structure will be presented in a forthcoming paper. 

V. ACKNOWLEDGEMENT 

We acknowledge support from ARO Grant No. W911NF-08-1-0317 and from EU through 
Grant No. PlTN-GA-2009-234970. 



1 M. I. Dyakonov and V. I. Perel, Phys. Lett. A 35, 459 (1971). 

2 J. E. Hirsch, Phys. Rev. Lett. 83, 1834 (1999). 

3 S. Zhang, Phys. Rev. Lett. 85, 393 (2000). 

^ S. Murakami, N. Nagaosa, and S.-C. Zhang, Science 301, 1348 (2003). 

^ J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, and A. H. MacDonald, Phys. Rev. 

Lett. 92, 126603 (2004). 
^ H.-A. Engel, E. I. Rashba, and B. I. Halperin, in Handbook of Magnetism and Advanced Magnetic 

Materials, edited by H. Kronmiiller and S. Parkin (Wiley, Chichester, UK, 2007), vol. V, pp. 

2858-2877. 
'^ R. Winkler, in Handbook of Magnetism and Advanced Magnetic Meterials, edited by 

H. Kronmiiller and S. Parkin (Wiley, Chichester, UK, 2007), vol. V, pp. 2830-2843. 
8 D. Culcer and R. Winkler, Phys. Rev. B 76, 245322 (2007). 
^ D. Culcer and R. Winkler, Phys. Rev. Lett 99, 226601 (2007). 

10 D. Culcer, E. M. Hankiewicz, G. Vignale, and R. Winkler, Phys. Rev. B 81, 125332 (2010). 

11 W.-K. Tse, J. Fabian, I. Zutic, and S. Das Sarma, Phys. Rev. B 72, 241303 (2005). 

12 A. Mal'shukov, L. Wang, C. Chu, and K. Chao, Phys. Rev. Lett. 95, 146601 (2005). 



20 



13 V. M. Galitski, A. A. Burkov, and S. Das Sarma, Phys. Rev. B 74, 115331 (2006). 

14 T. Tanaka and H. Kontani, New Journal of Physics 11, 013023 (2009). 

15 E. M. Hankiewicz and G. Vignale, J. Phys. Cond. Matt. 21, 235202 (2009). 
1^ G. Vignale, J. Supercond. Nov. Magn. 23, 3 (2010). 

1^ L. Liu, O. J. Lee, T. J. Gudmundsen, D. C. Ralph, and R. A. Buhrman, arXiv:1110.6846 (2011). 
18 M. W. Wu, J. H. Jiang, and M. Q. Weng, Phys. Rep. 493, 61 (2010). 

1^ J. Fabian, A. Matos-Abiague, C. Ertler, P. Stano, and I. Zutic, Acta Physica Slovaca 57, 565 
(2007). 

20 D. D. Awschalom and M. E. Flatte, Nature Physics 3, 153 (2007). 

21 I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 323 (2004). 

22 Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. Awschalom, Science 306, 1910 (2004). 

23 V. Sih, R. C. Myers, Y. K. Kato, W. H. Lau, A. C. Gossard, and D. D. Awschalom, Nature 
Physics 1, 31 (2005). 

24 J. Wunderhch, B. Kaestner, J. Sinova, and T. Jungwirth, Phys. Rev. Lett. 94, 047204 (2005). 

25 N. P. Stern, S. Ghosh, G. Xiang, M. Zhu, N. Samarth, and D. D. Awschalom, Phys. Rev. Lett. 
97, 126603 (2006). 

2^ N. P. Stern, D. W. Steuerman, S. Mack, A. C. Gossard, and D. D. Awschalom, Nature Physics 
4, 843 (2008). 

27 L. Liu, R. A. Buhrman, and D. C. Ralph, arXiv: 1111. 3702 (2011). 

28 S. Valenzuela and M. Tinkham, Nature 442, 176 (2006). 

2^ T. Kimura, Y. Otani, T. Sato, S. Takahashi, and S. Maekawa, Phys. Rev. Lett. 98, 156601 

(2007). 
3*^ T. Seki, Y. Hasegawa, S. Mitani, S. Takahashi, H. Imamura, S. Maekawa, J. Nitta, and 

K. Takanashi, Nat Mater 7, 125 (2008). 

31 P. Nozieres and C. Lewiner, J. Phys. (Paris) 34, 901 (1973). 

32 Y. A. Bychkov and E. I. Rashba, J. Phys. C 17, 6039 (1984). 

33 E. G. Mishchenko, A. V. Shytov, and B. I. Halperin, Phys. Rev. Lett. 93, 226602 (2004). 

34 R. Raimondi and P. Schwab, Phys. Rev. B 71, 033311 (2005). 

35 A. Khaetskh, Phys. Rev. Lett. 96, 056602 (2006). 

36 E. M. Hankiewicz and G. Vignale, Phys. Rev. Lett. 100, 026602 (2008). 

37 R. Raimondi and P. Schwab, Europhys. Lett. 87, 37008 (2009). 

21 



3* R. Raimondi, P. Schwab, C. Gorini, and G. Vignale, Ann. Phys. (Berlin) 524, 153 (2012). 

39 A. Ohtomo and H. Y. Hwang, Nature 427, 423 (2004). 

""^ S. Thiel, G. Hammerl, A. Schmehl, C. W. Schneider, and J. Mannhart, Science 313, 1942 

(2006). 
"^1 M. Huijben, G. Rijnders, D. H. A. Blank, S. Bals, S. V. Aert, J. Verbeeck, G. V. Tendeloo, 

A. Brinknian, and H. Hilgenkamp, Nature Materials 5, 556 (2006). 
"^2 E. Dagotto, Science 318, 1076 (2007). 
^3 A. D. Caviglia, M. Gabay, S. Gariglio, N. Reyren, C. Cancellieri, and J.-M. Triscone, Phys. 

Rev. Lett. 104, 126803 (2010). 
"^^ Z. S. Popovic and S. Satpathy, Phys. Rev. Lett. 94, 176805 (2005). 
'^^ H. Chen, A. Kolpak, and S. Ismail-Beigi, Phys. Rev. B 82, 085430 (2010). 
'^^ O. V. Dimitrova, Phys. Rev. B 71, 245327 (2005). 



22 



